Deterministic Preparation of a Tunable Few-Fermion System 
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Systems consisting of few interacting fermions are the building blocks of matter, with atoms and 
nuclei being the most prominent examples. We have created a few-body quantum system with 
complete control over its quantum state using ultracold fermionic atoms in an optical dipole trap. 
Ground-state systems consisting of 1 to 10 particles were prepared with fidelities of ~ 90%. We can 
tune the interparticle interactions to arbitrary values using a Feshbach resonance and observed the 
interaction-induced energy shift for a pair of repulsively interacting atoms. This work is expected 
to enable quantum simulation of strongly correlated few-body systems. 



The exploration of naturally occurring few-body 
quantum systems such as atoms and nuclei has been ex- 
tremely successful, largely because they could be pre- 
pared in well defined quantum states. Because these sys- 
tems have limited tunability, researchers created quan- 
tum dots — "artificial atoms" — in which properties such 
as particle number, interaction strength and confining 
potential can be tuned (7, 2). However, quantum dots are 
generally strongly coupled to their environment which 
hindered the deterministic preparation of well-defined 
quantum states. 

In contrast, ultracold gases provide tunable systems in 
a highly isolated environment (3, 4). They have been 
proposed as a tool for quantum simulation (5, 6) which 
has been realized experimentally for various many-body 
systems (7-10). Achieving quantum simulation of few- 
body systems is more challenging because it requires 
complete control over all degrees of freedom: the parti- 
cle number, the internal and motional states of the par- 
ticles, and the strength of the inter-particle interactions. 
One possible approach to this goal is using a Mott in- 
sulator state of atoms in an optical lattice as a start- 
ing point. In this way, systems with up to four bosons 
per lattice site have been prepared in their ground state 



(77, 72). Recently, single lattice sites have been ad- 
dressed individually (13). In single isolated trapping ge- 
ometries, researchers could suppress atom number fluc- 
tuations by loading bosonic atoms into small-volume 
optical dipole traps (14-18). However, these experi- 
ments were not able to gain control over the system's 
quantum state. 

We prepare few-body systems consisting of 1 to 10 
fermionic atoms in a well-defined quantum state making 
use of Pauli's principle, which states that each single- 
particle state cannot be occupied by more than one iden- 
tical fermion. Therefore, the occupation probability 
of the lowest energy states approaches unity for a de- 
generate Fermi gas, and we can control the number of 
particles by controlling the number of available single- 
particle states. We realize this by deforming the confin- 
ing potential such that quantum states above a well de- 
fined energy become unbound. This approach requires 
a highly degenerate Fermi gas in a trap whose depth can 
be controlled with a precision much higher than the sep- 
aration of its energy levels. 

To fulfill these requirements, we use a small volume op- 
tical dipole trap with large level spacing. This micro- 
trap is created by the focus of a single laser beam (fig. 
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Figure 1: Experimental setup. Systems with up to 10 
fermions are prepared with 6 Li atoms in a micrometer- 
sized optical dipole trap created by the focus of a single 
laser beam. The number of atoms in the samples is de- 
tected with single atom resolution by transferring them 
into a compressed magneto-optical-trap (MOT) and col- 
lecting their fluorescence on a CCD camera. A Fesh- 
bach resonance allows one to tune the interaction be- 
tween the particles with a magnetic offset field. 



1) with a waist of wq < 1.8 /mi and measured radial 
and axial trapping frequencies (cj r , uj a ) = 2tt x (14.0 ± 
0.1, 1.487 ± 0.010) kHz (79). We load the microtrap 
from a reservoir of cold atoms. The reservoir consists of 
a two-component mixture of 6 Li atoms in the two low- 
est energy Zeeman substates \F = 1/2, mp = +1/2) 
and \F = 1/2, m F = -1/2) (labeled state |1) and 
1 2)) in a large volume optical dipole trap. The reser- 
voir has a degeneracy of T/Tp « 0.5 (79), where Tp 
is the Fermi temperature. We superimpose the micro- 
trap with the reservoir and transfer about 600 atoms into 
the microtrap. After removal of the reservoir, the de- 
generacy of the system is determined by Tp ^3 /iK 
in the microtrap and the temperature T < 250 nK of 
the reservoir (20). Assuming thermal equilibrium be- 
tween the microtrap and the reservoir this corresponds 
to T/Tp < 0.08. According to Fermi-Dirac statistics 
this yields an occupation probability for the lowest state 
exceeding 0.9999 which is large enough not to constrain 
our preparation scheme. 

To spill the excess atoms from the microtrap we add a 
linear potential in the axial direction by applying a mag- 
netic field gradient. To obtain the same potential for 
both components we apply the gradient at a large mag- 
netic offset field where the difference in the magnetic 
moments of states |1) and |2) is negligible. A particu- 
lar magnetic field is then chosen so that the interaction 



strength between atoms in the different states vanishes 
because of a nearby Feshbach resonance (fig. S2). By 
varying the depth of the microtrap and the strength of 
the magnetic field gradient we can control the number 
of bound states in the potential (fig. 2A). If fewer than 
10 bound states remain, the system is essentially one- 
dimensional because of the approximate 1:10 aspect ra- 
tio of the trap; consequently, each energy level is occu- 
pied by one atom per spin state. During the spilling pro- 
cess we adiabatically tilt the potential, wait to let atoms 
in unbound states escape and then ramp the potential 
back up (79). 

To probe the prepared systems, it is necessary to mea- 
sure the number of atoms in the microtrap with single 
atom resolution and near unity fidelity. We achieve this 
by releasing the atoms from the microtrap, recapturing 
them in a compressed magneto-optical trap and then 
recording their fluorescence with a charge-coupled de- 
vice (CCD) camera (fig. 1). With this technique we can 
count the total number of atoms in the magneto-optical 
trap with a fidelity exceeding 99% for 1 to 10 atoms (fig. 
SI). 

Figure 2 shows the mean atom number and its variance 
as a function of the minimum microtrap depth during the 
spilling process. The atom number shows a step-like de- 
pendence on the trap depth with plateaus for even atom 
numbers. These plateaus appear at trap depths where 
the potential barrier for atoms in the uppermost level 
becomes so low that these atoms leave the trap on a 
timescale that is shorter than the duration of the spilling 
process. A simple estimation (27) shows that the life- 
time of this state can be up to three orders of magni- 
tude shorter than the lifetime of the lower states. When 
an appropriate trap depth is chosen the fluctuations in 
the atom number are as low as var/(N) = a 2 / (TV) = 
0.017 for eight atoms, corresponding to a suppression of 
1 8 dB compared to a system obeying Poissonian statis- 
tics. We can then calculate an upper bound for the de- 
gree of degeneracy in the microtrap of T/Tp < 0.19 
by assuming that all fluctuations result from holes in 
the Fermi distribution; this provides a complementary 
method to probe the degeneracy of the lowest energy 
states of an ultracold Fermi gas which is conceptually 
related to recent studies of antibunching in degenerate 
Fermi gases (22, 23). 

To estimate the probability of finding the system in its 
ground state after the spilling process we bin the mea- 
sured fluorescence signal into a histogram (fig. 3B). For 
the preparation of systems consisting of two fermions 
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Figure 2: (A) The spilling process. Starting from a de- 
generate two-component Fermi gas of about 600 atoms 
in the microtrap, we create few-particle samples by adia- 
batically deforming the potential to spill atoms in higher 
levels. After the potential has been restored, the system 
is in a well defined few-particle quantum state. (B) Con- 
trolling the number of quantum states. When the trap 
depth is reduced, the mean atom number decreases in 
steps of two because each energy level in the trap is oc- 
cupied with one atom per spin state. Each data point is 
the average of ~ 190 measurements with a as the stan- 
dard deviation and var = a 2 as the variance (shown on 
the right) (C). For even atom numbers, the number fluc- 
tuations are strongly suppressed. For eight atoms, we 
achieve a suppression of 18dB of var/(N) compared 
to a system obeying the Poissonian statistics. 



we obtain a fidelity of 96(1)%. The error is the statistical 
error calculated by assuming that the occurrence of sam- 
ples with undesired atom number follows a Poissonian 
distribution. From combinatorial considerations (79) we 
deduce that only a negligible fraction of the prepared 
two-particle systems are not in the ground state before 
we ramp the potential back up at the end of the prepara- 
tion process. To check whether we create excitations in 
the system by ramping up the potential, we perform the 
spilling process a second time. After the second spilling 
process we measure a fidelity of 92(2)% for preparing 
two atoms. This yields an upper bound of 6(2)% for the 
excitation probability during the potential ramps (79). 
If we assume the same excitation probability for ramp- 



ing up and down we get an estimated fidelity of 93(2)% 
to prepare the system in its ground state after ramping 
the potential back up after the first spilling process. For 
eight atoms we find a ground state preparation proba- 
bility of 84(2)%. By varying the time between the two 
spilling processes we found the 1 / e-lifetime of the pre- 
pared two-particle system in its ground state to be ~ 60 s 
which shows the high degree of isolation from the envi- 
ronment. 

To realize configurations with an arbitrary imbalance 
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Figure 3: (A) Fidelity of preparing systems in the 
ground state. To determine how many of the prepared 
few particle systems are in their ground state, we repeat 
the spilling process. This removes atoms in higher levels 
but leaves the ground state unchanged. (B) Histograms 
after the first and second spilling process for the prepa- 
ration of (i) two atoms and (ii) eight atoms. The num- 
bers above the peaks give the relative occurrences of the 
counts within the corresponding peaks. The fidelity af- 
ter the second spilling process (right) remains almost 
unchanged, indicating that the ground state is prepared 
with high fidelity. 

in the number of atoms in state |1) and |2) we prepare 
balanced systems and perform a second spilling process 
that only removes atoms in state |1). We do this by 
changing the value of the magnetic offset field to 40 G 
where atoms in state |2) have negligible magnetic mo- 
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ment and are therefore unaffected by the magnetic field 
gradient (fig. S2). Using this technique we have cre- 
ated imbalanced systems with fidelities similar to those 
in the balanced case (19). 

Precise control over the trapping potential is not only 
essential to prepare few-body samples, it is also an ef- 
fective tool to probe strongly interacting systems. We 
use it to explore one of the simplest non-trivial few-body 
systems: two non-identical fermions with repulsive in- 
teractions in the ground state of a one-dimensional har- 
monic trap. We first prepare a non-interacting pair of 
atoms in states |1) and |2) in the ground state of the 
trap. Then we lower the potential barrier such that 
the two atoms slowly escape the trap on a timescale of 
r = 630 ± 120 ms which we measure by recording the 
decrease in the mean atom number as a function of hold 
time (fig. 4B). 

To repeat the measurement for two atoms with repulsive 
interaction, we tune the scattering length a to a large 
positive value, a = 4100 ao, where ao is the Bohr ra- 
dius, using a Feshbach resonance. In this case we ob- 
serve a much higher initial loss rate followed by a slow 
decay. This fast initial loss can be explained by the 
energy shift U of the ground state due to repulsive in- 
teractions which effectively decreases the height of the 
potential barrier (fig. 4A). Because the system is quasi- 
one-dimensional for the lowest 10 axial states, one has 
to consider the radial confinement (24) for the calcula- 
tion of U (25). Given the trap parameters one expects 
a shift U on the order of half the axial level spacing 
per particle. After one of the atoms has left the trap, 
the interaction energy drops to zero, leaving the remain- 
ing atom in the unperturbed ground state of the poten- 
tial. Within our measurement accuracy, we measure an 
equal probability for the remaining atom to be in state 
|1) and 1 2). By developing theoretical models for these 
interaction-induced dynamics one can use this method 
to quantitatively study strongly interacting systems. 
The system we created is well suited for quantum sim- 
ulation with fully controlled few-body systems. For at- 
tractive interactions, it can be used to study BCS-like 
pairing in finite systems which is a model used for the 
description of nuclei (26). Splitting a trap containing 
a repulsively interacting pair of atoms into a double 
well creates entangled pairs of neutral atoms, which can 
be used for quantum information processing (27, 28). 
By transferring the prepared ground-state samples into 
a periodic potential (29), our system can be used to 
overcome one of the current major challenges of study- 
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Figure 4: (A) By tuning the barrier height, we can 
prepare systems where the atoms leave the trap with a 
well defined rate. Interactions cause an energy shift U 
and thereby change the effective height of the potential 
barrier. (B) The barrier is chosen such that two non- 
interacting atoms (a = 0) escape from the trap on a 
timescale of r = 630 ± 120 ms (black squares). With 
repulsive interactions (a = 4100 ao) we observe a much 
faster escape of one atom. After this atom has left the 
trap, the loss of the second atom becomes suppressed 
(red circles). Each data point is the average of ~ 190 
measurements. The error bars show the standard error 
of the mean. 



ing quantum many-body physics with ultracold atoms: 
preparing systems with sufficiently low entropy to ex- 
plore phenomena such as quantum magnetism. 
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Methods 

Preparation of the few-fermion samples. 

We prepare a reservoir of cold atoms by evaporative 
cooling of a spin mixture of states |1) and |2) at a 
magnetic field of B = 300 G as described in previous 
work (30). The reservoir consists of approximately 
2 • 10 4 atoms per spin state in a crossed beam optical 
dipole trap (A = 1070 nm, waist wo — 40 /im, crossing 
angle 15°) with trapping frequencies uj r « 2tt x 370 Hz 
and uji « 2tt x 34 Hz. According to a Gaussian fit 
to the in- situ density distribution the temperature in 
the reservoir is T < 250 nK, which corresponds to 
T/T F w 0.5. To load the microtrap we turn it on within 
100 ms and wait for 20 ms to let the sample thermalize. 
Then we ramp the magnetic field to B = 523 G, where 
the scattering length crosses zero to switch off the 
inter-particle interaction and subsequently remove the 
reservoir. 



Properties of the microtrap. 

The microtrap is created by the focus of a single laser 
beam (A = 1064 nm). It is focused with a two-lens ob- 
jective with a numerical aperture of 0.36. The beam 
waist in the focus is calculated from the measured trap 
frequencies to be wo ~ 1.8 /im. As the trap contains 
only few radial levels, the actual potential deviates from 
the harmonic approximation. To determine the energy 
levels in the trap for the lowest levels we prepare sys- 
tems of two non-interacting atoms in the lowest axial 
level of the trap. Then we modulate the microtrap posi- 
tion by shaking a mirror in the optical path of the micro- 
trap setup with frequency uj for a defined period. If uj 
equals to the difference in frequency of the first and sec- 
ond level AE12 / H, atoms are excited to the second level. 
We find AE 12 ,axiai/h = 2tt x (1, 487±0, 010) kHz and 
AE^radiai/h = 2tt x (14, ± 0, 1) kHz. 

Spilling scheme. 

To spill the atoms from the microtrap we first ramp up 
a magnetic field gradient of B' = 18.9(2) Gauss/cm 
within 150 ms such that approximately the ten lowest 
states remain bound in the trap. The maximum slope 
of the gradient is smaller than 1.5 Gauss/cm/ms which 
is slow compared to the frap frequency and therefore 



the sample is not heated during the ramp. The fine tun- 
ing of the potential is done 20 ms after the gradient has 
reached its final value by changing the depth of the op- 
tical trap. The spilling process consists of a linear ramp 
over 8 ms of the optical potential starting from 100% 
of the trap depth, corresponding to a laser power of 
P= 265(27) /iW, to 63-84% of the original depth (fig. 
2B), waiting for a hold time of 25 ms to let the atoms 
escape from the trap and linearly ramping the potential 
back to its original value within 8 ms. The potential of 
the tilted microtrap in the z-direction can be written as 

V(z) = V opt (z) + V rnag (z) 

= 1 '°< 1 -(TTW^)) ) - B '" 3 

where Vo /kb ~ 3 fjK is the maximum depth at the cen- 
ter of the optical trap, z r = ttwq/X is the Rayleigh 
range and /i is the magnetic moment of the atoms. 

High fidelity atom number detection. 

We detect the number of atoms in the prepared samples 
by recapturing them in a magneto-optical trap (MOT) 
and collecting their fluorescence signal (31). The mag- 
netic field gradient has a strength of 250 G/cm, the diam- 
eter of the MOT beams is ~ 4 mm and their frequency 
is red detuned from the resonance by twice the natural 
linewidth of the transition. While we cannot determine 
the recapture efficiency of the MOT directly, it cannot be 
lower than the highest measured preparation fidelity per 
atom of 98(1)%. To record the fluorescence from the 
MOT we image it onto a CCD camera with an imaging 
system with numerical aperture of 0.17, capturing about 
1% of the emitted photons. During the 0.5 s exposure 
time of the CCD one atom scatters about 1.9 x 10 6 pho- 
tons. Considering the numerical aperture and the quan- 
tum efficiency of the imaging system roughly 1 x 10 4 
photons per atom are detected. The exposure time is 
much shorter than the 1 / e-lifetime of 250 s of the atoms 
in the MOT measured for 8 atoms. This lifetime is long 
enough that neither light-induced collisions in the MOT 
nor collisions with background gas atoms limit our de- 
tection fidelity. To deduce the atom number from the 
fluorescence signal we bin all data from each series of 
measurements into one histogram. These histograms 
show distinct peaks, each corresponding to an integer 
number of atoms in the MOT. From the spacing of the 
peaks we extract the calibration factor for the mean flu- 
orescence per atom. Because of fluctuations of the in- 
tensity of the MOT beams and the detuning of the MOT 
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the fluorescence signal drifts on a few percent level on 
a timescale of several minutes. To compensate for this 
drift, we rescale the fluorescence signal of each mea- 
surement by a factor which is obtained by taking the av- 
erage fluorescence per atom of the ten previous and fol- 
lowing measurements. To obtain this average for rescal- 
ing we only consider data with fluorescence signals that 
are close to a peak in the histogram, i.e. maximum la 
distance. Then, the rescaled atom numbers are binned 
into a histogram (fig. SI) and Gaussians are fitted to the 
peaks. We find a separation of the peak centers of ~ 6a. 
The data points within 2a of a peak center are binned to 
integer values which represent the number of atoms in 
the prepared sample. The 5% of measurements outside 
of the two-sigma width of the peaks are rejected. This 
is possible because the measurements are uncorrelated 
and the atom number detection is independent from the 
preparation. 

Estimating the ground state preparation fi- 
delity. 

To estimate the fidelity for preparing two atoms in the 
ground state of the trap from the histogram shown in 
figure 4B we only consider the lowest levels of the trap. 
This is based on the assumption that atoms in higher lev- 
els have only negligible probability to remain trapped 
after the spilling process. 

From the number of prepared samples containing one 
or three atoms we can deduce upper bounds for the 
probability to find an atom missing in the lowest level 
or an atom remaining in the second level of 2% each. 
As the atoms are non-interacting during the preparation 
process we can assume these probabilities to apply to 
each atom individually. A system containing two atoms 
which are not in the ground state requires both a hole in 
the lowest level and an atom on the second level. This 
is suppressed by a factor of (0.02) 2 = 4 • 10 -4 . From 
this we conclude that only a negligible fraction of the 
observed two-atom samples were not prepared in their 
ground state. 

However this does not exclude the possibility of excit- 
ing the system when closing the trap at the end of the 
spilling process. Therefore we perform the spilling pro- 
cess a second time, which removes atoms in higher lev- 
els of the trap. We find that the second spilling process 
reduces the number of samples containing two atoms 
from 96(1)% to 92(2)%. If we assume that the 2% 
probability of preparing one atom after the first spilling 



process is due to states beeing non-occupied before the 
spilling process and account for the fact that almost all 
samples containing three atoms will have two atoms af- 
ter the second spilling process, we would expect 98(2)% 
of the samples to contain two and 2(2)% of the samples 
to contain one atom. This leads us to the conclusion that 
there is a 6(2)% probability to create an excitation while 
ramping the barrier up and back down. 
Following the same consideration for a system of eight 
atoms we find the probability that a sample of eight 
atoms was not in the ground state to be 4-1 -3 . From 
spilling twice we also find an upper bound of 6(2)% for 
the number of excitations during the ramps. 

Imbalanced samples. 

We create imbalanced samples with a different number 
of atoms in states |1) and |2) in a two-step process. 
We first prepare a balanced system and subsequently 
remove additional atoms in state |1) by performing a 
second spilling process at a magnetic field of 40 Gauss 
where the magnetic moment of state |2) vanishes (fig. 
S2). To estimate the preparation fidelity we consider 
the two components separately: The fidelity of the ma- 
jority component (|2)) is only given by the first spilling 
process since it is unaffected by the second spilling pro- 
cess. The uncertainty in the minority component (| 1)) is 
solely determined by the second spilling process. Thus 
the overall fidelity for the creation of an imbalanced sys- 
tem is similar to the balanced case. 
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Figure S2: Magnetic field dependence of the Zeeman 
sublevels and the scattering length. (A) Energy of the 
Zeeman sublevels of 6 Li in the F=l/2 electronic ground 
state. The inset shows the magnetic field region below 
100 G where the magnetic moment of state |2) crosses 
zero. (B) Two-body scattering length between atoms in 
states |1) and |2) as a function of the magnetic field. For 
evaporative cooling of the reservoir and for the transfer 
of atoms into the microtrap the offset field is tuned to 
300 G. The spilling process is performed at a magnetic 
field of 523 G where the interaction strength vanishes 
because of the zero crossing of the scattering length. For 
the observation of the interaction-induced dynamics we 
tune to the magnetic field 760 G to obtain a pair of re- 
pulsively interacting atoms. 



Figure SI: Histogram of the rescaled fluorescence 
signals from the measurements shown in figure 2. The 
black curves are Gaussian fits to the rescaled atom num- 
ber. The 2-atom (8-atom) peak center is separated from 
its adjacent peak centers by 7a (5.7a). This is large 
enough to clearly distinguish the fluorescence signal 
of different atom numbers. We bin fluorescence data 
within a 2 a- width of the peaks (red bars) to integer atom 
numbers; the counts outside the 2a widths(grey bars) 
are rejected. 
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